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COMPENSATION OF DISTRIBUTED DELAYS IN 
INTEGRATED COMMUNICATION AND CONTROL SYSTEMS 


ABSTRACT 

Advances in aircraft and spacecraft technology demand high-speed and 
reliable communications between the individual components and 
subsystems for decision-making and control. This can be accomplished 
by Integrated Communication and Control Systems (/CCS - ) which use 
asynchronous time-division-multiplexed networks. Unfortunately, 
these networks can introduce randomly varying distributed delays. 
This paper presents the concept, analysis, implementation, and 
verification of a method for compensating delays that are 
distributed between the sensor (s) , controller, and actuator (s) 
within a control loop. With the objective of mitigating the 
detrimental effects of these network- induced delays, a predictor- 
controller algorithm has been formulated and analyzed. 

Robustness of the delay compensation algorithm is investigated 
relative to parametric uncertainties in plant modeling. The delay 
compensator has been experimentally verified on an IEEE 802.4 
network testbed for velocity control of a DC servomotor. Dynamic 
performance of the delay compensator has also been examined by 
combined discrete-event and continuous-time simulation of the flight 
control system of an advanced aircraft, that uses the Society of 
Automotive Engineers (SAE) linear token passing bus for data 
communications. The paper is concluded with several areas of future 
research in the evolving field of ICCS . 
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1. INTRODUCTION 


Large-scale dynamical systems, such as advanced aircraft and 
spacecraft, require reduction of direct human intervention in 
control and decision-making processes and its replacement by 
hierarchical levels of automatic control [1] as much as practicable. 
In this respect, much effort is being focused on integration of 
expert systems, i.e., rule-based procedures, with model-based 
algorithms within the control loops. Such a multi-loop, intelligent 
control system requires extensive interactions between its disparate 
and spatially distributed subsystems and components. Such 
interactions, especially when the devices are not collocated, can be 
effectively carried out via an asynchronous time-division- 
multiplexed data communication network [2], This concept of Integrated 
Communication and Control Systems ( ICCS ) has already been adopted in 
aeronautics and astronautics [3], computer-integrated manufacturing 
(C/A/) [4], and is being actively pursued for integrated control of 
chemical processes [5] and future generation automobiles [6]. A 
general concept of how distributed delays are induced by an ICCS 
network is illustrated in Figure 1-1. 

The tasks of decision making, e.g. , fault detection, isolation and 
reconfiguration (FDIR) , at the system level are simplified within 
the common network environment of ICCS. The reason is that 
monitoring of status and dynamic behavior of the integrated control 
system becomes a relatively less difficult task due to availability 
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Figure 1-1. Example of an Integrated Communication and Control System 

























of the information that traverses throughout the communication 
network. This is of paramount importance in strategic processes like 
advanced aircraft and spacecraft. For example, a contingency plan 
may be implemented by reconfiguring the system software to 
compensate for failures. An additional advantage of using a network 
in FDIR is the convenient use of analytical redundancy (i.e., 
analytically deduced measurement (s) of a plant variable) to 
supplement the sensory information as well as for identifying common 
mode failures [7,8,9]. The above concept is also applicable to 
factory automation [10] in computer-integrated manufacturing (CIM) 
and integrated vehicle control systems in future generation 
automobiles. 

The network traffic in ICCS is heterogeneous in the sense that both 
re al~time (e.g., flight critical) and some of the non-real-time 
(e.g., avionics) functions are likely to share a common data bus. 
This would generate aperiodic and possibly random traffic in the 
network. The sensor and control signals in ICCS are subjected to 
randomly varying distributed delays although the average bus load is 
kept well below the saturation level [11-13]. These delays occur in 
addition to the sampling time and data processing delays. 
Furthermore, the ICCS could be subjected to recurrent loss of data 
due to noise corruption in the communication medium and malfunctions 
of the network protocol. Therefore, an observable (reachable) system 
that assumes availability of sensor (control) data at consecutive 
samples may become unobservable (unreachable) due to recurrent loss 
of data. If the plant is unstable in the open loop as it could be 



for highly maneuverable supersonic aircraft, then closed loop 
stability requires the system to be detectable and stabilizable. 
This is ensured by complete observability and reachability [14]. 
Recurrent loss of data could render the system undetectable and/or 
unstabilizable. Of importantance are the conditions for which this 
phenomenon could lead to loss of observability (reachability) 
[15,16] and therefore be a source of potential instability. This 
subject is discussed in depth in Appendix A. 


There are two parallel but complementary approaches to circumvent 

the problem of network-induced delays and loss of data: 

Approach #1 : Selection of a reliable communication protocol that will prioritize real-time 
traffic (e.g., sensor and control data for fast loops, and interrupt signals for event 
synchronization) by preempting the non-real-time traffic (e.g., aircraft systems operation and 
management data). An example is the fiber-optic version of the Society of Automotive 
Engineers (SAE) token bus protocol [17] that allows four levels of priority at a transmission 
speed of 100 Mbps and a high level of noise immunity. An appropriate protocol and a 
carefully designed network would certainly reduce the induced delays and the probability of 
data loss but their detrimental effects cannot be cjnpletely eliminated especially if the dynamic 
performance requirements of the control system are stringent or the bus utilization is not 
exceedingly low. ■ 

Approach #2 : Compensation of network-induced delays by: (i) modification of the existing 
control systems (that are designed for conventional digital control systems) by updating their 
parameters and incorporating additional observers/predictors, or (ii) development of a new 
control structure. This approach assumes that the distribution of the network traffic and the 
specifications of the access protocol are given, i.e., the statistical characteristics of the 
network-induced delays and the probability of data loss are assumed to be known from the 
point of view of ICCS design. ■ 


1.1. Benefits of the Reported Research 

Although significant efforts have been expended for improving the 
speed, reliability and other performance of computer networks, and 
ample research papers have been published in the areas of modeling 
and simulation of computer communication protocols, significance of 
network-induced delays and data loss relative to stability of 
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feedback control systems has not apparently received much 
consideration. Therefore, the approach #2 has been adopted in the 
reported research to deal with ICCS assuming that the network is 
well-designed using the approach #1. ICCS design via a combination 
of these two approaches will yield the following advantages: 

o Augmentation of operational flexibility and system reliability without increasing the network 
and computer hardware redundancy; 

o Enhancement of bus utilization, thereby, allowing incorporation of advanced control algorithms 
and Al-based decision support systems within the existing network and computer system; 

o Robustness of the integrated control system relative to noise, transient traffic, and malfunctions 
of the network, and uncertainties and disruptions in plant operation; 

o Reconfiguration of the control structure and self-repair (in the event of several unattended 
failures) via software migration over the network. 

1.2. Summary of Previous Research in ICCS 

The potential problem of network- induced delays in ICCS has been 
addressed in a sequence of three papers [11-13] under the NASA LeRC 
Grant No. NAG 3-823. In the first paper [11], a finite-dimensional 
model of ICCS is presented by taking into consideration the effects 
of network-induced delays. A necessary and sufficient condition for 
system stability was established for the special case of 

periodically varying (non-random) delays. Attention was focused on 
the control loops with identical sampling rates for the sensor and 
controller with a time skew A between the respective sampling 

instants. Although this model is more suitable for analysis of 

dynamic performance than other infinite-dimensional models [11], its 
application to stability analysis for non-periodically and randomly 
varying delays is not straight-forward. The second paper [12] 

discusses the issues related to ICCS design with an emphasis on how 
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A g can significantly contribute to network- induced delays and 

distort the control command sequence due to vacant sampling. An 

alternative approach, by which vacant sampling can be avoided, is to 

deliberately assign non-identical sampling periods T and T to the 

s c 

sensor and controller, respectively. The third paper [13] 

specifically addresses finite-dimensional modeling of control 

systems that have non-identical sensor and controller sampling rates 

and are subjected to randomly varying distributed delays. Criteria 

for selection of the sampling ratio T /T were presented for this 

s c 

special class of multi-rate sampling. 

1.3. Contributions of the Technical Report 

This report presents the concept, algorithm, and experimental and 
simulation results of a delay compensation algorithm for ICCS with 
identical sensor and controller sampling rates. The proposed delay 
compensator is intended to alleviate the detrimental effects of 
network-induced delays by using a multi-step predictor. The key idea 
in the compensator design is to monitor the data when it is 
generated and to keep track of the delay associated with it. With 
this knowledge, an algorithm keeps the delay constant as seen by the 
controller. Therefore, the closed loop control system model is 
constrained to be finite-dimensional, linear, time-invariant 
provided that the plant, observer, and controller are linear time- 
invariant. Therefore, this delay-compensated system model is more 
suitable for stability analysis of ICCS than the time-varying model 
reported in our previous publication [11]. 
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The major assumption in the above delay compensation algorithm is 
that the (randomly varying) network-induced delays are bounded. This 
assumption is justified in view of the fact that unbounded delays 
would render the feedback loop open. Using a specified confidence 
interval, an upper bound can be assigned to each of these randomly 
varying delays. The number of predicted steps in the compensator is 
then determined from these bounds, and the compensated system may 
use the original control law that was designed for the system 
without being subjected to the induced delays. The contributions of 
this report are: 

o Formulation of the concept of a delay compensation strategy for ICCS\ 

o Establishment of the separation principle for the multi-step predictor and state-variable 
feedback controller; 

o Analysis of some of the robustness properties of the delay compensator relative to plant 
modeling uncertainties; and 

o Implementation of the delay compensation strategy, and its testing by: 

- Experimentation with a servomotor control system at a network testbed that employs the 

IEEE 802.4 token bus as the medium access control protocol; 

- Simulation of the flight control system of an advanced aircraft using an SAE token bus 

network. 

1.4. Organization of the Technical Report 

The report is organized in seven sections and three appendices. The 
concept of using an observer for delay compensation and 
contributions of other researchers in this field are succinctly 
presented in Section 2. The delay compensation algorithm is derived 
in a closed loop representation in Section 3. Robustness analysis 
relative to modeling uncertainties is presented in Section 4 for a 
two-step delay compensator. A description of the experimental 



facility, details of the experiments, and interpretation of the 
experimental results are presented in Section 5. A functional 
description of the simulation program and a model of the simulated 
flight control system are presented in Section 6 along with a 
discussion and interpretation of the simulation results. The report 
is concluded in Section 7 with recommendations for future research 
to deal with the basic problems of stability analysis and robust 
design of the delay-compensated system. 

Appendix A presents the concept of extended observability in the 
event of recurrently random data loss, followed by pertinent 
results. Appendix B provides the supporting materials for proof of 
the delay compensation algorithm in Section 3. Stability of the 
delay-compensated system in presence of plant modeling uncertainties 
is discussed in Appendix C. 
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2. OBSERVER FOR ESTIMATION OF DELAYED STATES 


The sensor, controller, and actuator are assumed to have identical 
sampling rates. Furthermore, we assume the sensor and actuator to be 
synchronized, i.e., the time skew between their sampling instants is 
maintained at zero. Under these conditions, the plant and controller 
dynamics of the delayed system at the sample time k are modeled as: 

x k+1 = A + B -plant dynamics with delayed control (2-1) 

y k = C x k -sensor measurement (2-2) 

w k = y)c-A(k) -delayed output (2-3) 

u k = ^ W k^ -control function (2-4) 

where xe/? n , u eR s and ye./? 111 , and the matrices A, B, and C are of 
compatible dimensions; the finite non-negative integers T(k) and 
A(k) represent the numbers of delayed samples in control inputs and 
measurements, respectively. The control law u k is a linear function 
of the history W k := ( w j c / w ]c-i'”*} the delayed measurements. The 
plant in (2-1) and (2-2) is reachable and observable. The objective 
is to construct the control function £ such that the effects of the 
delays on the control system are mitigated. 

Remark 2J : It has been stated in Section 1 that the (randomly varying) network-induced delays are 
assumed to be bounded. Therefore, for the purpose of delay compensator design, T and A may be 
treated as upper bounds (in units of sampling periods) for the controller-actuator and sensor- 
controller delays, respectively. In contrast, the delays in the ICCS model reported in our previous 
work [11,12] are deterministically or randomly varying. ■ 
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Several investigators have addressed the problems of delay 
compensation in closed loop control systems. An intuitive approach 
[18] is to augment the system model to include delayed variables as 
additional states. Unfortunately, this renders some of the states 
uncontrollable even when the original system is completely 
controllable [19,20]. For the case of delayed control inputs, 
Pyndick [21] proposed a predictor for the optimal state trajectory 
based on past control inputs. Zahr and Slivinsky [22] considered the 
problem of controlling a computer-controlled system with measurement 
and computational delays. It was pointed out that the delays in 
multivariable systems may result in: (i) an increase in the 
magnitudes of the transients, and poor response during the inter- 
sampling time; (ii) loss of decoupling between individual SISO 
control loops although decoupling may be restored for a stable 
process at the steady state; and (iii) a possible decrease in 
stability margin. Their algorithm was verified by simulation but the 
use of an observer to estimate the unavailable states was not 
discussed. 

A significant amount of research work has been reported for observer 
and controller design [20,23,24] for the case of inherent constant 
delays that occur within the process to be controlled. In contrast 
to the system, under consideration in (2-1) to (2-4) where the 
sensor and control data are delayed, such processes are described as 
follows : 

= A x ( t ) + D x(t-h) + G u ( t) 

= c x(t) 


dx(t)/dt 
y (t) 
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(2-5) 

(2-6) 


where h is a constant. By setting G=0 and D=Br in (2-5) (where r is 
the gain matrix) , it reduces to a delayed state-variable-feedback 
system. 

The reported literature in delay compensation does not apparently 
address the problem of distributed delays in both the input and 
output variables, which is the case with ICCS . A possible approach 
[22] for compensation of constant delays that affect the input or 
output variables of a system is to predict the current output. 
However, if a state-space approach is used for predicting the 
output, the plant state variables must be obtained first. Moreover, 
the measurements might be corrupted with noise. In view of the 
above, we propose to use an observer for estimation of the delayed 
states and then to predict the current state using the state 
transition matrix. However, this multi-step observer would impose 
additional dynamics which generate additional phase lag in the 
control system. The impact of this phase lag on the system dynamic 
performance should not be significant because the controller is 
expected to be designed with a sufficiently large stability margin 
to allow for disturbances and uncertainties. 

The observer parameters can be optimally selected to attenuate the 
high-frequency noise in the measurements. If the measurement noise 
is considerable, then the proposed algorithm could be extended to 
include a stochastic filter [25] instead of a deterministic 
observer. Furthermore, a reduced order observer or a functional 
observer [26] can be used to reduce the observation lag. 



3. ALGORITHM OF THE OBSERVER-BASED DELAY COMPENSATOR 


The algorithm for compensation of distributed delays is built upon 
the concept of linear state-variable feedback control and multi-step 
prediction where no loss of sensor or control data is taken into 
account. The algorithm is presented below as a proposition. 


Proposition 3.1 : Given the following predictor-controller scheme: 


x k+i = Ax k + Bl V y k = Cx k Plant Model; (3-1) 
Z k+1 | 1 = Az k | 1 + Bu k + L k ( y k “ Cz k | Observer Model; (3-2) 
z k+l | r = Az k|r-1 +Bu k r-Step Predictor for r>2 ; (3-3) 
u k = r ] c z } c |p Predictive Control for a fixed p>2 ; (3-4) 


where z^| r : = estimate of x^ 

e k x k - z k|i is the 
the gain matrices (L^) 

the initial conditions 


given ( Y k _ r , Yjc-r-i » - ) ? (3-5) 

estimation error; (3-6) 

and (T^) are a priori determined, and 
and (u^, i=0, — ,p-l} are given. 


Then, the closed loop system equation can be expressed as 


> 


( 

\ 


r \ 

X k+l 


(A+Br k ) 

- Br k 3 k 


x k 

e k-p+l 


0 



e k-p 

J 


k. 

) 


^ J 


where the plant model is assumed to be exact and 



3 k :=< 


p p-1 p-i-1 

n(A-L k _ p+j _ lC ) + E [A^L^.C n(A-Vp + j-l C )] if P^ 2 

j=l i=l j=l 

if 0<p<2 


Proof ; The proof of this proposition is supported by lemma B.l in 
Appendix B. Subtracting (3-2) from (3-1) yields 


x 


k+1 


or 


z k+l | 1 
e k+l 


A(x k " z k|l ) V^k " z k|l ) 
(A - L k C) e k 


(3-7) 


Substituting (3-4) in (3-1) and using (B-l.l) from lemma B.l yields 

p ^ ^ 

*k + i - Ax k + B V z k|i -t 0 AiL k-i-i Ce k-i-i^ < 3 - 8 ) 

Adding and subtracting Br k x k and using (3-6) in (3-8) 

p^2 

x k +1 - < A+Br k) x k + Br k' e k -£ 0 AiL k-i-i ce k-i-i> < 3 - 9 ) 

Using (3-7) e k+g for some integer q>l can be expressed as 

e k+q “ tn i (A-L k+j _ 1 c )]e k (3-10) 

Replacing e k and in terms of e k _ p , i.e., the use of (3-10) in 

(3-9) yields 


x k+l ( A+BI V x k 
P 


kkk-p 

P- 2 i. 


where 3]c = [RA-L^.^C) + Yl A \-i-i C ? fljA-L^.^C) ] 
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The proof for p>2 follows by combining (3-7) and (3-11) . If p=0 
then 3 =1 since the delay is zero and the ordinary separation 

K 

principle is applicable. For p=l the ordinary separation principle 
still applies since only a first order prediction is used, i.e. u^ = 
r k z k|l* In other words / an ordinary linear state feedback controller 
with observer drives the plant using a first order prediction of the 
states; thus, a standard observer is naturally suited for 
compensation of a constant delay of one time step. Thus, if p=l then 

3 k =I * ■ 

Remark LL - If L.=L and r^r, i- e > constant observer and controller gains, then Proposition 3.1 
determines stability of the compensated system due to separation of the controller and the observer. 
That is, the eigenvalues of the closed loop delayed system are the same as the combined eigenvalues 
of the two matrices (A-LC) and (A+Br). ■ 

The delay compensator is implemented on the basis of a priori known 
upper bounds of the distributed delays that are induced by the ICCS 
network. These bounds are expressed as integer multiples of the 
controller s:aifipling period T (which is also the sensor sampling 
period) . The existence of these upper bounds is justified from the 
point of view of control system design because the control system 
will be open-loop otherwise. Figure 3-1 shows implementation of the 
delay compensation scheme in Proposition 3.1. The resulting control 
system is modeled as follows. 
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Figure 3-1. Schematic Diagram of the p-Step Delay Compensation Algorithm 














-delayed sensor data 


w 


k 


z k-p+l | 1 
z k+i | r 


u 


k 


y k-p 

(A-LC)z^,_p| 1 + Bu^ + Lw^ -observer dynamics 

Az k | r _ 1 + Bu^ for p>r>2 -p-step predictor 

rz, , -control law 

k|p 


(1-3) 

(1-4) 

(1-5) 

( 1 - 6 ) 


where xeR n , ueR s and y eR m , the matrices A, B, C, L, and r are as 
defined earlier. The number of predicted steps, p, in the 

compensator is determined, in general, as the sum of the upper 
bounds (relative to a specified confidence interval) of the 

controller-actuator and sensor-controller delays. 


Remark 3.2 : Loss of sensor and control data is not accounted for in the derivation of the delay 
compensation algorithm in Proposition 3.1. However, if the network recurrently loses data as it is 
pointed out in Section 1, the plant may be rendered unobservable or unreachable. Under these 
circumstances, the notion of observability (reachability) can be extended as explained in detail in 
[15,16] and Appendix A. If the plant satisfies the condition of extended observability (reachability), 
then the delay compensator may still function in a possibly degraded mode. Performance analysis of 
the delay compensator under recurrent loss of data is a subject of future research. ■ 


Remark 33 : The number, p, of predicted steps could be obtained as the sum of the specified bounds 
of the distributed delays, e.g., A(k) and T(k). If the joint statistics of A and T are known, p could 
be calculated more precisely. ■ 

Remark 3.4 : Although the sensor and actuator are assumed to be synchronized, the controller is not 
required to be synchronized with the sensor or actuator. The time skew between the controller and 
sensor/actuator sampling instants is absorbed within A and T. ■ 
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Remark 3.5: If the sensor and actuator are not collocated, then they can be synchronized by 
transmitting high priority interrupt signals via the network medium or by additional wiring [12], 
Nevertheless this will increase the system reliability requirements. ■ 
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4. MODELING UNCERTAINTY OF A TWO-STEP PREDICTOR/CONTROLLER 


In this section we investigate the effects of modeling uncertainty 
upon performance of the predictor-controller algorithm for the case 
when the total delay to be compensated is equal to two sampling 
periods. Let the plant model described in (3-1) have the true 
dynamical characteristics given as: 


k+1 


= A, 


+ B, 


'u. 


c, "x. 
k k 


(4-1) 


In view of the above the modeling uncertainties are defined as: 


6A, 


:= A k "-A; 


8 B, 


= V- B; 


<SC, 


:= C k "-C 


(4-2) 


Remark 4_A \ The postulated model of the plant is time-invariant whereas the uncertainties may be 
time-varying. ■ 

The controller structure, given by (3-1) -(3-6) with p=2 in (3-4) and 
steady-state controller gain r and observer gain and L, are: 


k 

= rz k | 2 

(4-3) 

k+l | 1 

" Az k|l + Bu k + L(y k ‘ Cz k|l> 

(4-4) 

k+1 1 2 

" Az k|l + Bu k 

(4-5) 

k 

= x k ” z k| 1 

(4-6) 


We intend to obtain a closed loop, augmented state variable 
representation in such a way that the modeled system dynamics and 
the modeling errors are decoupled to the maximum possible extent. 
From now onwards, for the sake of brevity, the subscript k from the 
time-varying matrices in (4-2) are omitted. 
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Substituting (4-1) and (4-4) into (4-6) yields 


k ~k-l + B " u k-1 [Az k-l| l +Bu k-l +L(C " x k-l Cz k-l|l )] 

k-1 “ * z k-l|l + 5Bu k-l 


= A"x 
= tf"x 

where := A"-LC" ; ^ := A-LC 


Adding and subtracting tyx^^ in the right hand side of (4-7) 

e k - * e k-i + ex k-i + <Bu k-i 
where © := 4'" -4' = 5A-L5C 

Substituting (4-3) in (4-8) yields 
e k 

Substituting (4-5) in (4-4) yields 

z k+l | 1 = z k+l|2 + Lc " x k " Lcz k | 1 


= »e k _ 1 + + «Brz k _ ll2 


'k| 2 *k | 1 LC " x k-l + LCz k-l|l 

= X k - < X k' Z k|l) - LC < X k-:T Z k-l|l> ' L(C"-C)x k _ 1 
Substituting (3-6) and (4-2) in the above equation yields 

z k | 2 = x k “ e k LCe k-l “ L5Cx k-l 
Substituting (4-9) in (4-10) yields 

z k | 2 = x k ^ e k-l + 0x k-l + 5Brz k-l| 2-* " LCe k-l " L5Cx k-l 
{ k ' (' t+ LC)e k _ 1 - (©+L6C) x k _ 1 - «Brz k _ l(2 


0 


= x. 


= X k - Ae k-1 - 5Ax k-l - SBrz k-l|2 


Substituting for z k _^j 2 in the above equation by using (4-10) 

z k | 2 x k Ae k-1 “ fiAx k-l “ 5Br ^ x k-l -e k-l -LCe k-2 -L5Cx k-2 ^ 

= x R - (A-SBF) e k-1 -5Gx k _ 1 +5BrL5Cx k _ 2 +$BrLCe k _ 2 (4-11) 

where SG := (<SA+SBr) 

Substituting for z k-1 | 2 in (4-9) by using (4-11) yields 

e k = ^ e k-i +0x k-l +6Br[x k-l" e k-2 _5Gx k-2 +5BrLiCx k-3 +5BrLCe k-3 

= ^ e k-l +5 B2 e k-2 +a B3 e k-3 +F A A B A C X k-l +a Bl X k-2 +a BC X k-3 (4 


(4-7) 


(4-8) 


(4-9) 


- 10 ) 


] 

- 12 ) 


13 



f a a b a c 

:= (0+6Br) = 

(5A-L5C+5Br) 

5 bi 

: =-5Br5G 

= -SBr(5A+SBr) 

a BC 

:= (6Br) 2 L5C 


3 B2 

: =-5Br (A-SBr) 


5 B3 

:= (5Br) 2 LC 



Substituting (4-12) into (4-11) yields 

Z k|2 - X k “( A-5Br ) t^ e k-2 +a B2 e k-3 +5 B3 e k-4 +F A A B A C X k-2 + 

3 Bl X k-3 +3 BC X k-4l ' SGx k-l + 5BrL5Cx k -2 + 5BrLCe k-2 
Grouping similar terms 


Z k | 2 X k 5Gx k-l + a A A B A C X k-2 + a B4 X k-3 + a B5 X k-4 


where 


+(5 B6 A ^ )e k-2 + 5 B7 S k-3 + a B8 e k-4 


= (A-5Br)5Br(«A+5Br) 

= - (A-$Br) (JBr) 2 L5C 


= (A-5Br)SBr (A-5Br) 

= -(A-5Br) (5Br) 2 lc 
Substituting (4-13) and (4-3) in (4-1) yields 

X k + 1 = A ” X k + B " r t X k " 5GX k-l + a A A B A C X k-2 + a B4 X k-3 

+ a B5 X k-4 + ( ' 9 B6 -A ^ )e k-2 + 3 B7 e k-3 + ^B8 e k-4^ 
Rearranging and grouping similar terms 


5 A a B a C 

:= - (A-5Br) F 

5 B4 : = 

-(A-6Br)a B1 

5 B5 : = 

-(A-5BD3 bc 

3 B6 : = 

5Br^+5BrLC 

3 B7 : = 

- (A-5Br) dg 2 

a B8 

-(A-5Br)3 B3 


X k+1 Yx k + °A A B X k-l + °A A B A C X k-2 + °Bl X k-3 + °B2 X k-4 


+ QS k-2 + °B3 e k-3 + D B4 e k-4 


(4 


where Y := (A"+B"r) = (A+Br)+£G 

Q := B'T(a B6 -A^) = B'T (5BrA-A (A-LC) ) 

D A a B :=_B " r * G =-B"r (SA+6Br) 

°A A B A C := B " ra A A B A C = B ” r [^ Br )L5C - (A-iBr) (5A-L£C+5Br) ] 


-13) 


-14) 
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D ! 

D I 

D ! 

D r 


B1 

:= B " ra B4 

B2 

:= B»ra B5 

B3 

:= B"ra B7 

B4 

:= B " ra B8 


= B"r(A-5Br)6Br(«A+6Br) 
=-B"r ( A-6Br) (6Br) 2 L5C 
= B"r(A-5Br)<5Br(A-5Br) 
=-B'T(A-5Br) ( 5 Br) 2 LC 


The closed loop equation of the 2-step delayed system follows from 
(4-12) and (4-14) 


x k+l 


Y D 

x k 


I 0 

x k-l 


0 I 

x k-2 

= 

0 0 

x k-3 


0 0 

e k-l 


0 0 

e k-2 


0 0 

e k-3 


0 0 


'A A B °A A B A C °B1 °B2 Q 


°B3 °B4 


0 

0 

I 

0 


0 

0 

0 

I 


f a 

A A B A C B1 


0 

0 

0 

0 

a 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 


0 

0 


BC 
0 I 

0 0 


« a 


B2 B3 
0 0 

I 0 


k-1 


X 


k-2 


k-3 


k-4 


'k-2 


'k-3 


'k-4 


(4-15) 


J 


where the subscripts of elements are chosen to describe their 
relative dependence on the respective errors in formulating the 
plant model matrices. The above equation can be separated into the 
nominal and uncertain parts: 


x k+ i - (V + SV) x k 


(4-16) 


T T T T T T T T n T 

! k X k-1 x k-2 X k-3 x k— 4 e k-2 e k-3 e k-4 J ■ 


V is the nominal part which is constant as it contains the 
constant matrices A, B and C of the plant model in (3-1), and 

5V is the uncertainty part which is possibly time-varying 
because of modeling error matrices 5A, 6B and SC. 


Stability of the above system is discussed in Appendix C. 
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Remark 4.2 : The modeling errors SA, SB and SC constitute a four-fold increase over the order of the 
augmented system matrix in (4-15). However, by setting SB=0, the following simplifications result. 


5 B1 5 B2 3 B3 5 BC °B1 °B2 °B3 °B4 0 and Y ( A+ * A+Br ) » D a a B _ BF5A ' 
D A-B-C =_BrA(5A-L5C) ' Q=~BrA(A-LC) , F A ^ C = (6A-LSC) . 

Accordingly the system is reduced with SB=0 as follows 


k+1 


k-1 


'k-1 


A+SA+Br -BWA -BrA (SA-LSC) -BrA(A-LC) 
10 0 0 

0 10 0 

0 0 SA-LSC A-LC 


k-1 

C k-2 

e k-2 


(4-17) 


Remark 4J \ By separating the uncertainty part, the above equation can be expressed in the format 
of (4-16) with reduced-order state representation (due to the assumption SB=0) as follows: 


**k + i = + ^ A 


(4-18) 


where ^ = [x^ ^ T e k _ 2 V, and 


V = 


s 



■\ 


S’ 




A+Br 

0 

0 

-BrA (A-LC) 


SA 

-BrsA 

-BrA (SA-LSC) 

0 

i 

0 

0 

0 

' x 5V - 

0 

0 

0 

0 

0 

I 

0 

0 


0 

0 

0 

0 

0 

0 

0 

A-LC 


0 

0 

SA-LSC 

0 




j 





J 


Remark 4.4 : With 5B=0 the effects of SA and SC are to increase the order of the closed loop system 
by two. ■ 
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5. EXPERIMENTAL VERIFICATION AT A NETWORK TESTBED 


The delay compensation algorithm was verified at a network testbed 
which uses the 10 Mbps IEEE 802.4 linear token passing bus and IEEE 
802.2 as the link layer protocol in an ISO compatible network 
architecture (see Appendix B of [12]). Three terminal interface 
units, each of which are installed as a pair of cards in the 
backplane of a microcomputer, serve as the communication link in the 
network testbed. These host microcomputers can communicate with 
each other via the network medium or directly by using RS-232C 
communication ports (at a maximum data rate of 9.6 kbps). These 
microcomputers are PC-AT compatible under the Disk Operating System 
(DOS ) . They are designated as host #1, host #2, and host #3. 

The host #3 serves as the network manager. By use of the remaining 
terminal interface units, namely host #1 and host #2, a traffic load 
generator (TLG) has been designed to emulate the scenario of a large 
number of virtual stations with varying traffic. These virtual 
stations form a virtual ring in the sense that a virtual token hops 
around this ring. All odd-numbered virtual stations are emulated on 
host #1 and all even-numbered virtual stations on host #2. These two 
microcomputers communicate with each other via the network. While 
the real token travels back and forth between hosts #1 and #2, a 
virtual token hops around this virtual ring to emulate an actual 
logical ring consisting of a large number of terminals. Although 
this arrangement does not guarantee the transmission rate of the 
actual network, it suffices to emulate the randomly distributed 
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delays between the source and destination terminals. The TLG can 
emulate any number of virtual stations, and the message inter- 
arrival time and message length at each station can be arbitrarily 
selected using the random number generator (RNG) functions that are 
resident in hosts #1 and #2. In this way, the TLG offers flexibility 
of generating randomly varying network-induced delays with specified 
upper bounds . 

In the test apparatus, the plant and actuator are represented by a 
constant-field DC servomotor, and a tacho-generator serves as the 
sensor for measuring the motor speed. The sensor and actuator are 
connected to one of the virtual stations in host #1 which provides 
the hardware and software for the DASH-16 A/D and D/A conversion via 
direct memory access (DMA) . This odd-numbered virtual station is 
designated as #1 in the emulated logical ring. The controller is 
located in host #2 and can be designated as any even-numbered 
virtual station, i.e., the controller station can be arbitrarily 
located in different positions in the logical ring relative to the 
sensor/actuator station. Since the host #2 functions as an integral 
of the TLG, multitasking operations are needed for concurrent 
execution of the TLG and controller functions. To avoid this 
multitasking operation (which is not convenient under the DOS 
environment) , the host #3 (which is practically idle during normal 
operations in the capacity of a network manager) is directly 
connected (via RS— 232C) to host #2 such that host #3 serves as a 
coprocessor for host #2. The control algorithm is resident in host 
#3 which receives the sensor data (originally generated by the A/D 
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converter at host #1) from host #2 and transmits back the control 
command to host #2. This control command is eventually transmitted 
via the network medium to host #1 for being fed to the D/A converter 
as shown in Figure 5-1. 

Remark 5-7: The transmission delays of the sensor and control data via RS-232C between hosts #2 
and #3 can be viewed as part of the total data processing delay at the controller. ■ 


Remark 5-2 : The time skew between the sensor and controller can be set to any desired value 
(between 0 and T) via the RS-232C link between hosts #1 and #3 during the initialization of this 
test facility. ■ 

The servomotor serves as a continuous-time plant which communicates 
with its computer controller via the network. This represents a 
control loop within the ICCS. A relative measure of efficacy of the 
proposed delay-compensation algorithm is obtained by comparing the 
performance of a Proportional-Integral (PI) controller without any 
delay compensation against that of the same controller with delay 
compensation . 


The steady-state characteristics of the plant were experimentally 
determined as: 


r 



K (u 


ss 


+ 0 ) 


for u >0 
ss K 

for -0<u ss <0 

for u <-£ 
ss 


(5-1) 
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Op-Amp 



Figure 5-1.. Schematic Diagram of the Experimental Facility 











where y gs = measured steady-state angular velocity of the motor, 
u gs = constant input voltage, 

K = steady-state plant gain, and 0 = limit of dead band. 


In the linear region, the motor dynamics were represented by a first 
order model after appropriate bias compensation in the control 
input: 


x 


k+1 


a x k + b u k 


(5-2) 


y 


k 


c x, 
k 


(5-3) 


where a = exp(-T/r), T and r being the sampling period and the motor 
time constant, respectively. 


(j+l)T 

b = J exp(-( (j + l)T-t)/r) K dt/r = [ 1-exp (-T/r) ] K, and c=l. 
jT 


The plant parameters were identified from experimental data to be: 
t=1.7s; K = 3.33 (rad/s) /volt ; and 0=3.4 volts which is equivalent 
to 11.3 rad/s. The state variable feedback, in this first-order 
plant, yields a proportional controller where the gain, 7 , is 
obtained by pole placement such that the closed loop pole p=a-h^ . 
Similarly the observer gain t is chosen to place the observer pole 
at f=a -t. 


Since the plant model is not exact, the observer is expected to 
generate a steady-state error in the state estimate. Therefore, a 
weak integral action was used to compensate for this possible 
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steady-state error. The integral gain was set, by trial and error, 
between 0.09 and 0.15 volt/(rad/s) to provide zero steady-state 
error for a variety of induced delays. 

Two different experimental setups were considered. The first 
experimental setup was oriented towards determining the efficacy of 
the delay compensation algorithm when the total delay in the control 
loop is a constant with magnitude equal to an integer number of 
sampling periods. The second experiment emulated an ICCS environment 
by using randomly varying delays. The objective of the second 
experiment was to examine the robustness of the observer-based 
compensator due to data rejection and vacant sampling [11,12] 
coupled with plant modeling uncertainties and nonlinearities. 

5.1. Experimental Setup #1: Velocity Control with Constant Delays 

The objective of this experiment was to investigate the performance 
of the uncompensated and compensated control systems of the 
servomotor with deterministic delays such that a constant bound of p 
sampling periods could be strictly enforced. This was accomplished 
by setting: (i) the number of virtual stations to a minimum of two; 
(ii) network-induced delays to desired values by adjusting the 
emulated traffic; and (iii) the time skew between (co-located) 
sensor and actuator sampling instants to zero. This arrangement is 
equivalent to generating a constant delay by storing the incoming 
data from the A/D converter into a FIFO, push-pop buffer of size p. 
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Remark 5-4 : From Proposition 3.1 in Section 3, it follows that the observer should behave identically 
for p=0 and p=l for constant gains T and L provided that the plant model is exact. ■ 

Remark 5-5 : Since the plant under test is nonlinear, the dynamic characteristics of the control 
system are a function of the reference input. ■ 

Figure 5-2 compares the response of the compensated control system 
with that of the control system without the observer when the delay 
was zero, and a step change of 16 rad/s in the reference velocity 
was applied from 0 rad/s. The controller and observer poles were set 
at p^O.Ss 1 and f=0.15s" 1 , respectively, the integral gain of the PI 
controller to 0.15 volt/ (rad/s) , and the sampling period to 0.3s. 
The maximum input to the motor was restricted to 10 volts to avoid 
any potential damage; this corresponds to a steady state velocity of 
about 22 rad/s. 

As expected, the performances of the system with and without the 
compensator were comparable for zero delay. In this case, the 
observer-based controller suffered from the error in state 
estimation. The tracking properties of the observer can be improved 
by feeding in well-structured information about the plant 
nonlinearities. For instance, the input to the observer can be 
modified as follows: 



r 

u 

con 

- 0/K 

if 

u > P/K 

con y/ 


U ob = 1 

0 


if 

1 u 1 </3/K 

1 con 1 ~ H/ 

(5-4) 


u 

con 

+ 0/K 

if 

u <- 8/K 

con H/ 





29 




(oas/pej) A^tdoxsa jex n6u V 


30 


Figure 5-2. Dynamic Response of the Motor: 

Fixed Delay = 0, Reference Velocity = 16 rad/sec 


is the 


where u . is the modified input to the observer and u 

ob c con 

actual input to the DC motor. This empirical modification mimics the 
nonlinearities in the steady-state characteristics of the servomotor 
by adding a priori known information about the bias to the observer in 
the form of an input, and does not increase the order of the 
observer. Results, to be presented next, indicate that the steady 
state tracking errors were practically eliminated in this modified 
system. With the above modifications, the controller and observer 
poles were reset at p=0.1s _1 and f=0.0, respectively, the integral 
gain to 0.09 volt/ (rad/sec) , and the sampling period was retained at 
0.3s. 


Remade 5-< 5 : The above parameter settings were tested, by trial and error, to yield best responses of 
this nonlinear system for different delays. ■ 

Remark 5-7 : In general, setting the observer pole at zero, i.e. making the observer to be deadbeat, is 
not a recommended practice because of possible noise amplification. The integral action, enforced in 
this experiment, should act a low pass filter to attenuate the high frequency noise. The response of 
the compensated system was found to settle down slightly faster than that of the uncompensated 
system when the delay was zero in both cases [17]. A possible explanation for this phenomenon is 
that a lead action is provided by the observer with a pole at zero. ■ 

Experiments were conducted for different combinations of delays and 
reference velocities. Series of curves comparing the responses of 
the compensated and uncompensated system are recorded in [17]. We 
summarize the general information first and then present the 
selected results in detail. 
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The performance of the compensated system was always superior to 
that of uncompensated system for delays of T and 2T, i.e., one and 
two sampling periods. For delays of 3T and larger, the compensated 
system generally exhibited significant oscillations, and the 
uncompensated system performed worse. For example, with a delay of 
3T at the reference velocity of 17 rad/s, the compensated system 
response exhibited a limit cycle which was 3 times smaller in 
magnitude than that in the uncompensated system. The impact of 
reference input on the nonlinear system is presented next. 

Because of the deadband nonlinearity of the servomotor, a sudden 
jump in the system response was caused whenever the reference input 
was set close to «12 rad/s. As the reference input was increased, 
the dynamic responses of both compensated and uncompensated systems 
improved regardless of the magnitude of the delay. This is expected 
because as the reference speed is increased the motor excitation is 
drawn away from the on-off nonlinearity region. Moreover, the plant 
model is possibly better represented at higher velocities where the 
motor voltage versus speed characteristics are more closely follow a 
straight line pattern. As a result of a dead band in the motor 
characteristics, limit cycles frequently appeared in the 
uncompensated system response. 

Now we present selected results in Figures 5-3 and 5-4 to show 
comparisons of dynamic responses of the delay compensated system for 
delays of T and 2T with step changes of 15 rad/s and 13.5 rad/s, 
respectively, in the reference velocity from 0 rad/sec. The response 
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Figure Dynamic Response of the Motor: 

Fixed Delay = T, Reference Velocity = T5 rad/sec 
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Figure 5-4. Dynamic Response of the Motor: 

Fixed Delay = 2T, Reference Velocity =13.5 rad/sec 


of the compensated system is less oscillatory and evidently superior 
to that of the uncompensated system in Figure 5-3. The uncompensated 
system in Figure 5-4 suffered from sustained oscillations because 
the reference velocity of 13.5 rad/s is close to 12 rad/s when the 
servomotor tends to stall. In contrast, the compensated system 
brought the response to the reference point without any noticable 
undershoot and thus saved the system from entering into the deadband 
zone. 

Remark 5-fl : The predictive properties of the compensator apparently largely eliminate the problems 
of limit cycling for delays of T and 2T. Since the above observation relates to a first order plant 
under specific operating conditions, these experimental results do not provide a definite conclusion 
regarding robustness of the delay compensator relative to plant modeling uncertainties and 
nonlinearities in higher order plants. ■ 

5.2. Experimental Setup #2: Velocity Control with Random Delays 

The objective of this experiment was to investigate the performance 
and robustness of the proposed delay compensation algorithm when the 
communication network was subjected to random traffic generated from 
a large number of virtual stations. The control system is likely to 
be subjected to data rejection and vacant sampling at the 
controller's receiver buffer [11,12] if the compensator is not used. 
This problem can be apparently circumvented by an appropriate choice 
buffer size, namely p, at the controller and actuator buffer. 
However, if p is selected to be small, the randomly varying delays 
would exceed this bound and generate vacant samples with a larger 
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probability. On the other hand, if p is selected to be large, the 
impact of randomly varying delays would be reduced at the expense of 
potential instability due to plant modeling uncertainties and 
nonlinearities. (Note: If the plant is unstable in the open loop, 
then it follows from Figure 3-1 that the dynamic error of estimation 
will be worse with a larger p) . Selection of p becomes a design 
problem which calls for trade-off between mitigating the effects of 
uncertainties in network traffic and plant modeling. 


Setup #2 is similar to setup #1 with the following exceptions: 


Number of virtual stations = 30. 

Identical traffic distribution at all 28 virtual stations except the sensor and controller. The 
message arrival process is Poisson with an average inter-arrival time of 0.3s. These virtual 
stations were assumed to have identical message lengths which were adjusted to keep the 
offered traffic [6] at 0.5. 

In view of additional noise due to randomly varying delays, locations of controller and 
observer poles were moved back to p=0.5s 1 and f=0.15s , respectively, from the previous 
values of 0.1 and 0.0. The integral gain was reset to 0.15 volt/( rad/s) from 0.09 to compensate 
for potentially larger steady-state estimation errors. 

In general, the system dynamic response deteriorated under random 
traffic for both compensated and uncompensated systems [5]. A 
typical comparison is shown in Figure 5-5 where the response of the 
compensated system is seen to be superior to that of the 
uncompensated system. Since the traffic distribution is Poisson, the 
upper bounds of the network-induced delays can only be set with a 
certain confidence. Having the bounds of each of the two induced 
delays set to the sampling period T implies that p=2 . This causes 
some of the sensor and controller data to be lost, and yields a 
steady-state error possibly due to vacant sampling. If it was 
possible to increase T without appreciably affecting the system 
dynamics, upper bounds of the induced delays would have larger 
confidence intervals. Consequently, the probability of exceeding 
these bounds would be reduced resulting in improved performance of 
the compensated system. 

Remarjk 5-9 : Both compensated and uncompensated systems were less sensitive to data rejection and 
vacant sampling, when the controller pole was set closer to unity. On the other hand, bringing the 
controller pole close to unity results in a sluggish system response. Design of the compensated 
control system under randomly varying delays needs a trade-off between these two opposing 
requirements. ■ 
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Figure 5-5. Dynamic Response of the Motor: 

Random Delay = 2T , Reference Velocity = 18 rad/sec 



6. SIMULATION OF A FLIGHT CONTROL SYSTEM 


The delay compensation algorithm was verified by simulation of the 
longitudinal motion control of an advanced aircraft. The flight 
control system model and the combined continuous-time and discrete 
event simulation program were similar to those reported in [2). 
Simulation results were generated by augmenting the continuous-time 
part of the above model with the proposed delay compensation 
algorithm while the discrete-event model of the network was 
unchanged . 

Simulation results were generated under the following conditions: 


The network protocol is the optical version of SAE linear token passing bus with a 
transmission rate of 100 Mbps. 

The network has 30 stations that share the common medium. 

Station #1 operates as the sensor terminal, and station #2 as the controller terminal with its 
transmitter queue handling actuator commands and its receiver queue handling sensor data. 

Terminals #1 and #2 have a sampling period of 50ms and constant message lengths with the 
information part equal to 64 bits. 

Time skew (between sensor and actuator sampling instants) to zero. 

Random traffic in stations #3 to #30 was identically distributed. The message arrival process 
is Poisson with an average inter-arrival time of 50ms. These terminals are assumed to have 
identical message lengths which are adjusted to maintain the offered traffic [27] at any desired 
level. 


A schematic diagram of the flight control system is shown in Figure 
6-1. (A discrete version of the controller was implemented in the 
simulation.) The state-variable model of the plant, excluding the 
sensor data conversion factor of 180/rr, is described below. 
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PLANT 



Figure 6-1. Block Diagram of the Flight Control System 
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Plant Variables and Parameters: 


6 = Elevator command, i.e., input to the actuator (radian) 

cl 

$ e = Elevator deflection, i.e., actuator output (radian) 

W = Normal component of linear velocity 

at the center of mass (ft/s) 

Q = Pitch rate about the center of mass (radian/s) 

a = Angle of Attack (radian) 

A = Normal component of linear acceleration 
n at the sensor location (ft/s 2 ) 


The dimensional stability derivatives [2,28] for longitudinal motion 
dynamics were selected as: 


Z, 

de 

: = 

(dZ/dSJ/m 

= -202.28 ft/s 2 . 

Z 

q 

: = 

(dZ/3Q)/m 

= -16.837 ft/s, 

z 

w 

: = 

(dZ/dW ) /m 

= -3.1332 s' 1 . 

M de 

: = 

(dM/d6 e )/I y 

= -40.465 s' 2 . 

M q 

: = 

(3M/9Q)/Iy 

= -2.6864 s' 1 . 

M 

w 

: = 

(5M/aW)/I y 

= -0.01429 (s-ft)' 

M . 
wd 

: = 

(dM/dW)/I y 

- -0.00115 ft' 1 . 


where M is the pitch moment; 

Z is the normal component of the aerodynamic force; 
m is the lumped mass of the aircraft; and 
I is the moment of inertia about the pitching axis. 



Other constant parameters were: 


g = 32.2 ft/s 2 


(Acceleration due to gravity) 


l = 12.268 ft 
t = 0.05 s 


(Distance between the center of gravity of 
the airframe and the accelerometer) 
(Actuator time constant) 


U Q = 1005.3 ft/s (Reference flight speed) 


Longitudinal Motion Dynamics in the Continuous Time Domain : 


dx/dt = A x + B u; y = C x (6-1) 

where x = [* e W Q] T , u = 5 , y = [a A n Q] T , and 
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[M +M . (Z +Z ) 1 . 
L q wd v q w' J ' 



< M d S +M wd Z de> ' s 2 : = <V M wdV' 

- < Z de + ' S l )/g ' s 5 := <V' S 2 )/9 ' 


and 


S 6 := (Z g +' S 3>/«- 


The sensor and controller sampling period T was chosen to be 50ms in 
contrast to 10 ms in [2]. The rationale for selecting a larger value 
of T was to illustrate how the delay-compensated system would 
perform relative to the uncompensated system when the sampling rate 
may have to be reduced to allow for longer processing time for 
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executing complex control algorithms. This problem might be 
encountered in future generation hypersonic aircraft. 

The observer in the delay compensator was designed in a transformed 
state space where the three measured outputs a, A n , and q were made 
the state variables. All three poles of the discrete-time observer 
were located at 0.3s 1 . Then the compensator was constructed 
following the scheme in Figure 3-1 and the control law in Figure 6- 
1. In the simulation program, the conversion factor of 180/* as 
shown in Figure 6-1, was used after the estimated states were 
computed by the observer. 

The dynamic responses of the controller output, angle of attack, 
normal acceleration, and pitch rate were obtained for a unit step 
change in the reference signal under different offered traffic in 
the network. The compensated system response was always less 
oscillatory than that of the uncompensated system. From the point of 
view of aircraft control, the benefits of the less oscillatory 
response are better dynamic response, and reduced control effort and 
actuator wear. 

Figure 6-2 presents dynamic responses of pitch rate for the three 
cases of non-delayed (i.e., without network-induced delays), delay- 
compensated, and uncompensated delayed control systems when the 
offered traffic was set to 0.5. Apparently, the detrimental effects 
of network-induced delays are partially circumvented and the flight 
control system response is improved significantly with the use of 
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Figure 6-2. Dynamic Response of Pitch Rate 


the compensated control scheme. Further research is needed for 
robustness evaluation of the observer relative to plant modeling 
uncertainties beyond the analysis presented in Section 4. 

Remark 6^J_ : A common tendency in ICCS design is to increase the sensor 
sampling rate so as to avoid large delays introduced by the 
communication network on the tacit assumption that an arbitrarily 
large bandwidth is available in optical data communication networks. 
The above assumption may not hold if the network is shared by a 
large number of control and decision-making functions. Eventually, 
the increase of network traffic may result in large induced delays 
for some of the control loops even though the average bus load may 
remain well below the bandwidth capacity. Following a similar 
argument, the number of subscribers served by the network would be 
unduly restricted if the network- induced delays have to be bounded 
within a very conservative limit. Nevertheless, if appropriate 
models of the controlled processes are available then, as seen in 
the experimental and simulation results, the network utilization can 
be increased by using compensated control systems instead of 
increasing the sampling rates. The additional computational 
requirement for the compensator, depicted in Figure 3-1, is expected 
to be small in comparison to that of the complex algorithms of 
integrated control systems in the current and future generation 
advanced aircraft. Therefore, this concept of observer-based delay 
compensation will allow the controller computers to have less 
stringent processing power requirements. ■ 
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7. SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS FOR FUTURE RESEARCH 

A procedure for compensating the effects of distributed network- 
induced delays in Integrated Communication and Control Systems 
( ICCS ) has been proposed. The problem of analyzing system^ with 
time-varying and possibly stochastic delays could be circumvented by 
use of a deterministic observer which is designed to perform under 
certain restrictive but realistic assumptions. 

The proposed delay-compensation algorithm is based on a 
deterministic state estimator and a linear state-variable feedback 
control law. The deterministic observer can be replaced by a 
stochastic observer without any structural modifications of the 
delay compensation algorithm. However, if a feedforward-feedback 
control law is chosen instead of the state-variable-feedback control 
law, then the observer needs to be modified in the same way a 
conventional non-delayed control system should be. Under these 
circumstances the delay compensation algorithm would be accordingly 
changed. 

The separation principle of the classical Luenberger observer [26] 
holds true for the proposed delay compensator. 

The proposed delay-compensation algorithm is suitable for Integrated 
Communication and Control Systems (ICCS) in advanced aircraft, 
spacecraft, manufacturing automation, and chemical process 
applications. If the individual components of the ICCS communicate 
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with each other via a common communication medium, then the network 
can be designed such that the induced delays are bounded. In this 
way the detrimental effects of distributed delays induced by the 
network can be circumvented by the compensator. 

As the first step to robustness analysis, the impact of modeling 
uncertainties, i.e., errors in the plant model matrices, on the 
performance of the delay-compensated system has been investigated. 
Further analytical research is recommended following the established 
techniques for robustness and stability [29-32]. 

This report has addressed the concept, analysis, implementation, and 
verification of an algorithm for compensation of network-induced 
delays that occur in Integrated Communication and Control Systems 
(ICCS ) . As reported in Section 3, this algorithm has been developed 
on the basis of a deterministic state estimator, a multi-step 
predictor, and a linear state-variable feedback control. The 
proposed delay-compensation algorithm is suitable for ICCS in large- 
scale processes like advanced aircraft, spacecraft, autonomous 
manufacturing plants, and chemical plants. The communication 
network, that interconnects the individual subsystems and components 
of the ICCS, should be designed such that the induced delays are 
bounded relative to a specified confidence interval. In this way the 
detrimental effects of the distributed and varying network-induced 
delays can be (at least partially) circumvented by the compensator. 
The impact of modeling uncertainties, i.e., errors in the plant 
model matrices, on dynamic performance and stability of the closed 


46 



loop compensated system has been investigated for a two-step delayed 
system. Given the specifications of the control system, the above 
approach has a potential for establishing bounds on the modeling 
errors. 

Implementation of the delay-compensation algorithm, followed by 
experimental and simulation results have been presented. 
Experimentation with a DC-motor in a network testbed assembly has 
demonstrated that the observer-predictor functions in presence of 
plant modeling uncertainties, nonlinearities and disturbances, and 
measurement noise. Simulation of the flight control system of an 
advanced aircraft within a network environment also shows that the 
delay compensator is capable of improving the system dynamic 
performance. However, further analytical and experimental research 
beyond those reported in this report is needed to enhance robustness 
of this delay compensator. 

Possible areas for future research in the evolving field of 
Integrated Communication and Control Systems are innumerable. Some 
of the topics related to the research reported in this report are 
furnished below. 

Research ToqJc. #/ : Stochastic Filtering 

A stochastic filter can be designed to accommodate the effects of 
varying delays. The following formulation can be used to specify the 
filter requirements. Let the plant model be described as 
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x k + l = A k x k + B k u k + n k 


*k = c k x k + v 


where the plant and measurement noises, n^ and v^, are assumed to be 
white, mutually independent, and of covariances Q and R, 
respectively. The delayed sensor data w^, as input to the filter, is 
represented as: 


w k 5 0^k )Y k + 5 l^k^ y k-l + 5 2^ (p k^ y k-2 + *” + 5 p-l y k-p+l 


or, alternatively, as 


w k = WVk + s l (‘’k* c k x k-l + •" + s p-l (ip k )c k x k-p+l + v k- 
where <S^(j), i , j=0 , 1, 2 , ••♦, p-1 is the Kronocker delta defined as 

1 if i=j 

S- (j) = , and 

0 if i<j 

the network- induced parameter <? k is a stochastic chain defined as: 


( k, cj) : Kxf2 — ► { 0 , 1 , 2 , ••• , p—1 } , 


with K and fi representing the above finite set of non-negative 
integers and the sample space, respectively. 


The objective is to obtain the estimate, x^, of the state x^, which 
minimizes the performance index 



where is the collection of past measurements { w g, w i 

E { • | } is the conditional expectation given W^, and 
M is a positive semi-definite weighting matrix. 
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The above algorithm, formulated with p=2, can be applied to the 
problem of sequential state estimation of discrete processes with 
interrupted observations. Sawagari et al . [33] proposed a similar 

formulation but the input to the filter was formulated, with p=2 , 
as : 


W k 5 0 ((P k )C k x k + 


'k* 


In the case of missing data, a filter based on the proposed 
formulation may yield better results than that obtained by Sawagari 
et al . Particularly, the approach of Sawagari et al. may not be 
effective if x^ is under steady state conditions because, assuming 
that v^ is small, w^. would oscillate between 0 and C^x^. On the 
other hand, in the proposed formulation, w^ can vary only between 

c k-i**-i and c k x k- 


Research Topic #2 : Robust Compensator Design using LQG/LTR 

The performance robustness and stability of the control system would 
be affected by insertion of the delay compensator in the loop. 
Having designed the optimal state feedback gain r, the observer gain 
L in the compensator needs to be tuned so that the overall system 
performance is optimal under plant modeling uncertainties, process 
disturbances, and measurement noise. Linear Quadratic Gaussian (LQG) 
using Loop Transfer Recovery (LTR) [31,32] is a potential candidate 
for designing the delay compensator. 
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Research Tonic #3 : Optimization of Control and Network Parameters 


A cursory treatment on optimization of network- induced delays and 
associated parameters has been given in Section 4 of [12] for two 
special configurations of the control system: One is the case of 
identical sampling frequency of the sensor and controller within a 
feedback control loop; the other addresses the situation where the 
sensor sampling frequency is larger than that of the controller, 
which is a viable option for ICCS design. 

For identical sampling, the critical network parameter is A, the 
time skew between the instants of sensor and controller samplings, 
which is a very slowly varying parameter. Therefore, A should be 
periodically adjusted to maintain optimal performance of the delayed 
control system. The control law designed with the assumption of a 
fixed delay may not perform satisfactorily when the system is 
subjected to network-induced delays. Furthermore, as the 
characteristics of network traffic change, the induced delays would 
vary. Consequently, the controller parameters and the time skew A 
should be updated. Therefore, the objective is to derive an optimal 
control law for a system that is subjected to network— induced delays 
under varying network traffic. Although research addressing optimal 
control has been reported for systems with transportation lag 
[34,35], the simultaneous optimization of control and network 
parameters, namely the state-variable feedback matrix F and time 
skew A, has not been apparently attempted before. Future research 
may make use of the model derived in [11]. 
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For the alternative approach of non-identical sampling [12.13], the 
research problem is identification of an optimal sampling ratio e, 
which minimizes a specified performance index and guarantees 
stability of the closed loop control system under given network 
traffic statistics and a fixed controller sampling period. 

Research Topic # 4 ; Selection of the Number of Predictor steps 

The delay compensation algorithm, presented in this report, does not 
take into account the effects of recurrent loss of sensor and/or 
control data in the network. (See Remark 3.2 and Appendix A). The 
compensated system is expected to perform in a gracefully degraded 
mode if the system remains observable and reachable under recurrent 
loss of data. The number, p, of predicted steps may be required to 
be increased to accommodate this feature. Selection of p is dictated 
by uncertainties in both network traffic and plant modeling. 
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APPENDIX A 


EXTENDED OBSERVABILITY UNDER RECURRENT LOSS OF OUTPUT DATA 

An observable (reachable) system that assumes availability of sensor 
(control) data at consecutive samples may become unobservable 
(unreachable) due to recurrent loss of data in the computer network. 
If the control system is unstable in the open loop as it could be 
for highly maneuverable supersonic aircraft, then recurrent loss of 
data could render the system undetectable and/or unstabilizable . 

The concept of state estimation with recurrent loss of sensor data 
has been addressed by other investigators in different contexts. For 
example, Sawagari et al . [33] and Jaffer and Gupta [36] considered 
the problem of sequential estimation with interrupted observations 
within a stochastic setting. The sequential state estimation 
algorithms in both cases are developed using a Bayesian approach. To 
the best of the authors' knowledge, the problem of observability 

under persistent and random loss of data has not been studied 
before . 

This appendix introduces the concept of extended observability under 
recurrent loss of data where the state vector has to be 
reconstructed from an ensemble of sensor data at non— consecutive 
samples. Given that the computer network is designed to keep the 
probability of losing more than m data in every set of (i/+m) 
successive data less than an a priori prescribed bound p , the problem 
is to establish test criteria for observability (reachability) under 
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this condition, v is the observability (reachability) index, and m 
is known as a function of p and u. A fixed-size window of (iH-m) data 
from the available collection is considered for each observation, 
and sample numbers of the missing data are routinely recorded by the 
computer network protocol. This concept is different from that of 
extended observability with unknown inputs, for example, Basile and 
Marro [37], Emre and Silverman [38], Kudva et al. [39], Molinari 
[40], and Rappaport and Silverman [41]. 

The extended observability can be determined by testing the rank of 
every possible matrix associated with an augmented set of output 
data arising from each possible combination of data loss. Although 
this test is exhaustive, it is time-consuming, in general, and could 
lead to incorrect conclusions due to computational inaccuracy. An 
alternative approach for determining extended observability in the 
single-output case is presented. The relevance of the work reported 
in this appendix is summarized below: 

o A necessary and sufficient condition for extended observability which can be expressed via a 
recursive relation. 

o Necessary conditions (for the above) that are related to the characteristic polynomial of the 
state transition matrix in a discrete-time setting, or of the system matrix in a continuous-time 
setting. 

o A system-theoretic approach for having an insight into the problem of loss of observability. 

A.l. Extended Observability: Concepts and Test Criteria 

Let the plant be represented by a discrete-time, linear, time- 
invariant model in a deterministic setting at the sampling instant k 
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x ( k+ 1 ) = A x (k) + B u (k) ; z(k) = C x(k) 


(A-l) 


where the state vector xeR n , the input vector ueJ? r , the output 
vector z€f? p , and the constant matrices A, B, and C are of 

compatible dimensions. Furthermore, rank of C is p and the pair 
(C,A) is observable with observability index v. Then, 

i j_1 i 

x(k+j) = A J x (k) + Y. A B u(k+j-l-i) (A-2) 

i =0 

Defining y(k+j)= CA^xCk) , it follows from (A-l) and (A-2) that 

j_1 i 

y (k+ 3 ) = z (k+j ) - [ C A B u (k+j-l-i) (A-3) 

i =0 

The modified output vector sequence (y(k) ) can be used for state 
reconstruction in lieu of (z(k)} provided that the input sequence 
(u(k)) is available. (Observability with unknown inputs [37,41] is 
not addressed here.) Under normal circumstances, the state x can be 
reconstructed from v consecutive sets of outputs. However, in the 
event of loss of outputs, y(k) is not available at every k. The 
problem is to determine whether the state x(k) can be reconstructed 
by selecting any u output vectors from the collection Y m (k) of (i/+m) 
output vectors as defined below. (Note: Y m (k) may also be viewed as 
a (p(j/+m)xl) vector.) 

Y m (k) = [y(k) T | y (k+l) T |y(k+2) T | ... | y (k+v+m-1) T ] T 

where m is a fixed finite integer. 

Definition A. I : The system in (A-i) is said to be m-observable if x(k) can be reconstructed from any 
v distinct vectors in Y m (k). ■ 
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We now elucidate the concept of m-observability by second order, 
single output systems which are observable but not m-observable. 


Example AJ_ : Let A and c be defined as follows. 


Given A 


© 



0 

1 

fl 



and c = [1 0], the observability matrix 


implies that the pair (c,A) is observable. 


cA 


0 -1 


2 

On the other hand, cA = [-1 0] = -c. 

The system is not 1-observable and hence, cannot 
m>l. From a geometric point of view, A is a 90° 
implying that cA* +2 =-cA* for i= 0,1,2,3«*. ■ 


be m-observable, 
rotation matrix 


Example. AJ. '- Next we consider a 120° rotation matrix so that 


A 


- 1/2 -^ 3/2 

V 3/2 - 1/2 


c = [1 0 ]. 


This system is 1-observable but not m-observable, m>2 . ■ 

Next we proceed to determine the conditions for m-observability. The 
following relationship between Y m (k) and x(k) is derived from (A-l) , 
( A-2 ) and (A-3). 


Y m (k) = Qx(k) ( A-4 ) 

where Q = [© T | © m T ] T which is a (p(i^+m)xn) matrix, 

0 = [(CA°) T | (CA 1 ) 1, | | (CA i/ “ 1 ) T ] T , and 

© m = [(CaV | (CA^V | ... | (CA" +m - 1 ) T ] T 


55 



We now establish a simplified necessary and sufficient condition for 
m-observabil ity . 


Proposition AJ_ : The system described in (A-l) is xn-observable in the 

single output case (i.e., C is a lxn row vector c and y=n) iff all 

minors of © 0 1 are non-zero, 
m 


Proof : It suffices to show that an arbitrary set of n rows from Q© 
contains linearly independent vectors iff all minors of © m © 1 are 


non-zero. Consider the (n+m)xn matrix 



I n is the (nxn) identity matrix and V=© m © 


(A- 5) 


Let us choose any (n-k) rows from I and any k rows from V. If k=0, 
or k=n for m>n, this proposition automatically holds. Therefore, we 
consider l<k<n. Let U represent the (n-k) -dimensional subspace 
spanned by the (n-k) rows selected from 1^ and let W represent the 
subspace spanned by the k rows selected from V. First suppose that 
all minors of V are non-zero. Having all k-th order non-zero minors 
implies that any k rows are linearly independent and each linear 
combination of these rows must have fewer than k zero elements, 
i.e., more than (n-k) non-zero elements. Therefore, dimW=k and 
WnU={0}. Consequently, any n rows of Q© ^ are linearly independent. 
This establishes sufficiency. Next suppose that V contains a zero 
minor of order k. Then it is possible to construct a linear 
combination of these rows with at least k zero elements. Such a 
vector can be expressed as a linear combination of (n-k) rows of I n> 
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Therefore this collection of k+(n-k)=n rows of Q© 1 does not form a 
linearly independent set. This establishes necessity. ■ 


Next we relate certain properties of the matrix A in (A-l) to m- 
observability for the single-output case. 


Observation AJ_ : The matrix © m © ^ is completely determined by the 
coefficients of the characteristic polynomial of A. This can be seen 
by expressing © m © 1 as follows. 


-1 n T 

© m © x = [(cAV 


. ,n+l.T 
(cA ) 


(cA n+m - 1 ) T ] T ©- 1 


= [ (cA n © 1 ) T (cA n+1 © 1 ) T 


( cA n+m_ 1 ©~ 1 ) T ] T 


= [(c© 1 ©A n © 1 ) T (c© 1 ©A n+1 e“ 1 ) T 


(c©" 1 ©A n+in ~ 1 © 1 ) T ] T 


r / A n \ T / T 

= [ (cA ) (cA ) 


, .n+m-l.T,T 
(cA ) ] 


(A- 6) 


where A = ©A© -1 and c = c© -1 , 


This represents a similarity transformation into the standard 
observability canonical form [26] as follows. 


where A 


0 

0 


1 

0 


' a 0 " a i 


0 

1 


0 

-a. 


-a 


n-1 


Therefore (A-6) depends only on 
characteristic polynomial of A. ■ 


and c = [1 0 0 . . . 0] ( A-7 ) 


the coefficients of the 
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Observation A. 2 : Following (A-6) and (A-7) in Observation A.l, the first 

row of 0 G" 1 is [-a_ -a. ••• -a . ] . By Proposition A.l, a 

m '•0 1 n-l 

necessary condition for m-observability , m>l, is that all lxl minors 

of © e" 1 be non-zero. This implies that each coefficient of the 

m 

characteristic polynomial of A must be non-zero for m-observability, 
m>l . ■ 


Observation AJ : The characteristic polynomial of A in Example A.l is 

A 2 +I , i.e., a^=0. Therefore, the system is not m-observable Vm>l. 

But, the characteristic polynomial in Example A. 2 is A 2 +A+I, 
implying that the necessary condition for m-observability is 
satisfied. ■ 

Observation AA - The degree of m-observability is common to all 

observable pairs that share the same state transition matrix. That 
is, if 3 c" such that (c" ,A) is m-observable, then (c,A) is m- 

observable for each observable pair (c,A) . The reason for the above 
is that the output vector c is only required to establish the 

invertibility of 0, i.e., observability. On the other hand, given 

that 0" 1 exists, © m © 1 depends solely on the coefficients of the 

characteristic polynomial of the A-matrix. ■ 

The test procedure in Proposition A.l requires computation of the 

matrix V=0 0 1 . We now show how the coefficients of the mxn matrix V 
m 

can be recursively computed. From (A-6) and (A-7) it follows that 
the coefficients of V can be expressed as 

V(i,j) = V(i-l,j-l) - V(i-l,n) a j _ 1 l<i<m and l<j<n 
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by setting the initial conditions V(.,0)=0 and V(l,j)=-aj_ 1 . Then, 

n-1 

V(i, n) = V( i-1 , n-1) -V(i-l , n) a = -£ a.V(i-n+k,n) 

k=0 K 

Defining = V(i,n), the above equation can be expressed as 
n-1 


V i k 5 0 a k V i-n+k 


for i>n 


The range of V.^ in (A-8) can be extended 
V x = -a n _ 1 and defining V Q = 1 and = 0 


(A-8 ) 

by using the relationship 
for i<0 as follows. 


V. 

1 


k£ 0 akVi - n+k 


for i>l 


(A-9 ) 


Using (A-6) and (A-7) the i— th row of V can be expressed as 

[V(i,l) V ( i , 2 ) ... V(i,n)] = cA n 1+ '*' for i>l (A-10) 

Hence, from (A-10) and (A-ll) , it follows that cA n+i , for i>0, can 
be expressed as 


cA 


n+i 



V. 

l 


0 


n-1 


] 


0 


0 


where V . 

i 


< 




1 


k ?o akVi_n+k 


v i-l *• 

V. , , . 

l- (n-1) 

v i •• 

v. . 0 , 

l- (n-2 ) 


0 .. v i J 

for i < 0 
for i = 0 

for i > 0 


(A-ll) 


Extended observability is critical for design of networked control 
systems [11-13] especially if the plant under control is unstable in 
the open loop. This concept is also applicable to bad data 
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suppression [42] that may lead to random rejection of signals in 
feedback control systems. 

A natural extension of the reported work is to establish test 
conditions for extended observability in multiple-output systems. 
This requires the loss of data to be considered from two different 
perspectives: Case 1: The block of data containing the (pxl) output 
vector z(k) is unavailable; and Case 2: Only certain elements of z(k) 
are unavailable at the instant k while the remaining elements could 
be used for state reconstruction. In the latter case identity of the 
unavailable outputs may vary with k. 
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APPENDIX B 


Lemmas for Proposition 3.1 


This Appendix presents the four lemmas that are necessary to derive 
the delay compensation algorithm in Proposition 3.1 in Section 3. 


Lemma B.l: 


r-2 

: k|r = z k|l -X Q AlL k-i-l Ce k-i-l for 142 


(B-l.l) 


Proof : The proof of lemma B.l requires lemmas B.2, B.3, and B.4. We 
also introduce a definition that will be required to establish 
recursive relations in the lemmas. 


Definition BA : f k := , g R := z } 


G (k, j ) := 


W i 

£ A g k-i + l for 


otherwise 


F(k,j), : 


A* fk - i+1 

0 otherwise 


for k>j>l 


By use of Definition B.l and lemma B.2, it follows that 


z k I r = G (k-1 , r-1) - F (k-1 , r-2 ) 


(B-l. 2) 


Using the expression for G(.,.) in Definition B.l in conjunction 
with lemma B.3, i.e., the equation (B-3.1)), we have 


r-2 . 


r-2 . 


G(k-l,r-l) =.1^.1 


r-2 . 


r-2 . 


= f L Alf k-i 


(B-l. 3) 
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Using the expression for F (• , - ) in Definition B.l in conjuction with 
(B-l. 3) , it yields 

r-2 . 

G(k-1, r-1) = f k + F (k-1, r-2) - ^^..^Ce^.^ (B-l. 4) 

The proof follows by substituting (B-l. 4) and the expression for f k 
(from Definition B.l) in (B-l. 2). ■ 

Lemma B.2 : Following Definition B.l, 

z k+l I r = G (k,r-1) - F ( k > r " 2 ) (B-2.1) 

Proof: The identities G(k,l)=g k+1 and F(k,0)=0 are obtained from 
Definition B.l, respectively. Using (B-l. 2) in conjunction with 
these identities yields 

Z k+l|2 = G (k, 1 ) - F (k, 0) (B-2.2) 

Using (3-3) from Section 3 with r=3 and r=2 and subtracting 


z k+l | 3 z k+l | 2 A(z k|2 z k|l } 

Using lemma B.3 and (B-2.2) in (B-2.3) 

Z k+1 1 3 = G (k,l) ” F ( k '°) + A (9 k " f k ) 

= (G(k, 1) + Ag k ) - (F(k, 0) + Af k ) 

= G ( k, 2 ) - F(k, 1) 


(B-2.3) 


(B-2.4) 
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The method of induction is now used to complete the proof of 
Proposition 3.1. Using (3-3), as in (B-2.3), results in 


: k+l | r+1 z k+l | r A(z k|r z k|r-l ) 


(B-2.5) 


Substituting (B-2.1) in (B-2.5) 


z k+l | r+l -z k+l | r = A t G (*-l,r-l) - F(k-l,r-2) - G(k-l,r-2) 

+ F (k-l, r-3) ] 

= A[G (k-l , r-1) -G (k-l , r-2 ) ] -A[F (k-l , r-2 ) -F (k-l , r-3 ) ] 


Setting j=r-l in (B-4.1) of lemma B.4 and similarly j=r-2 in (B- 
4.2), and then substituting these results in the above equation, we 
obtain 


z k+l | r+l -z k+l | r = t G ( k ' r > “ G(k,r-1)] - [F(k,r-1) - F(k,r-2)] 

= [G (k, r) - F (k, r-1) ] - [G(k,r-1) - F(k,r-2)] (B-2.6) 

Using (B-2.1) in the right hand side of (B-2.6) yields 

z v+1 | r+1 = G(k,r) - F (k, r-1) ■ 


Lemma B.3 : 


z k | 2 z k | l ^k— l Ce k— l 


or equivalently, g k = f R - L k-i Ce k-l 


(B-3 . 1) 
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Prool : From (3-1), (3-2) and (3-6) in Section 3, it follows that 

z k|i = Az k-i|i + Bu k-i + VA-i < b - 3 - 2 > 

Also, using (3-3) , we have 

z k|2 = Az k-l|l + Bu k-1 • ( B - 3 - 3 ) 

The proof follows by substituting (B-3.2) in (B-3.3). ■ 

Lemma B.4 : 

G(k,j+1) - G (k, j ) = A[G (k-1 , j ) - G(k-1, j-1) ] (B-4.1) 

F(k, j+1) - F (k, j ) = A[F(k-l,j) - F(k-1, j-1) ] (B-4.2) 


Proof i 


j i j_1 i 

G(k, 3 + 1)-G(k,j) = £ A g. _ - £ aV f 


i=0 K +1 i=0 


D" 1 i D-2 i 

G (k-1 , ] ) -G (k-1 , j -1) = £ A x g . - £ A x g. . 

i=0 K 1 i=0 K_1 


(B-4.3) 

(B-4.4) 


(B-4.1) is obtained by substituting (B-4.4) into (B-4.3). The proof 
of (B-4.2) follows a similar procedure. ■ 
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APPENDIX C 


Stability of the Delay-Compensated System 


This appendix discusses stability of the delay-compensated system in 
presence of plant modeling errors, as described by (4-15) and (4-16) 
in Section 4. The objective is to outline an approach for 
determining a sufficiency condition for stability of the time- 
varying systems in (4-16) and (4-18) . 

In this sequel we select a norm of the nominal matix V such that 
this norm is close to the spectral radius of V. Thus, if the 
uncertainty is zero, i.e., SV=0, a measure of stability margin is 
provided by this norm that would be close to that of an eigenvalue 
test. This provides a necessary and sufficient condition for 
stability of linear time-invariant systems. 

Given a matrix Ve/? nxn and any e>0 there exists a norm IMI*/ which is a 
function of both V and e such that p (V) < || V || (V) +e where p(V) is the 

spectral radius of V, i.e., the magnitude of the largest eigenvalue 
[43]. The norm is obtained as follows. 

Transform V into its Jordan form J such that V = PJP’ 1 where P is an 
appropiate transformation matrix. Construct a diagonal matrix D 
with the i th diagonal element equal to e 1-1 , 0<i<n, and let Z = PD. 
The desired norm is || V || ^ = ||Z ^VZfl , that is, the infinity norm of 
transformation D a JD. Therefore, 
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(C-1) 


II V I * =max £ e j ‘ 1 |J i1 l 
i j=l 

The above equation circumvents the problem of matrix inversion for 
large n and small e. 

It follows that if p(V)<l we can select a small e and use the 
following inequality to obtain a test for closed loop stability 

» x k+ l«* = I (V+SV(k))X k ||* < ||V+fiV(k)|* ||X k ||* (C-2) 

Since ||V+SV(k)||* = | Z'Vz+Z'^V (k) Z | and Z'Vz is a constant matrix, 
the computation burden is on calculating Z _1 5V(k)Z only. 

The selection of the norm ||*||* is an attempt to reduce the amount of 
conservatism that is inherent in norm tests such as those using the 
natural norm. Indeed, it is possible for a matrix with spectral 
radius smaller than one to have an Euclidian norm, or any other 
standard norm, greater than one. In addition, the set 

E* = {n: || V+n || ^ < 1} (C-3) 

constitutes a region in R nxn such that if 6V(k)eE A for all k, then 

the control system described in (4-1) -(4-4) is stable. However, it 

should be noted that that Oce^^ does not necessarily imply that 

E(V,e ) is contained in E(V,e ) or vice versa. 

2 1 
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